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Resumen 

En este trabajo, se presenta un modelo numerico para el analisis de paneles pianos de 
mamposterfa. El comportamiento mecanico de la mamposterfa se Simula mediante un modelo 
continuo elastoplastico danado para materiales homogeneos, que considera no solo la fisuracion, 
sino tambien la naturaleza ortotropa del panel. El analisis se realiza mediante un programa de 
elementos finitos 2D, con no-linealidad material. 

Teniendo en cuenta que la mamposterfa tiene una estructura periodica, las propiedades 
mecanicas de la misma, se obtienen del analisis numerico previo de una celda basica 
representativa del conjunto, que se discretiza con una malla densa distinguiendo cada material. 
Esto permite simular el comportamiento de la mamposterfa a partir de los materiales 
componentes, ladrillos y mortero, y de la disposicion geometrica que ellos ocupan. 

El trabajo se completa con ejemplos de aplicacion en paneles de mamposterfa simple y 
reforzada, que prueban la validez del modelo propuesto. 



Abstract 

A numerical model for the analysis of plane masonry panels is presented in this paper. The 
mechanical behavior of masonry is simulated through a continuum elastoplastic damage model 
for homogeneous materials that considers not only cracking but also the orthotropic nature of the 
panel. The analysis is made with a material non-linear 2D finite element program. 
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Taking into account the periodic structure of masonry, the mechanical properties are obtained 
from a previous analysis of a basic cell. A fine mesh is used for this cell that is considered to be 
representative of the whole panel. 

In this way, departing from the component materials and their arrangement, the behaviour of 
masonry can be simulated. 

The paper is completed with application examples in plain and reinforced masonry panels that 
validate the proposed model. 



1 INTRODUCTION 



La mamposterfa es hoy, particularmente en el mundo en desarrollo, el principal material 
estructural para resolver las crecientes demandas urbanas. Existe en todo el mundo, gran cantidad 
de edilicios de mamposterfa o de hormigon armado combinados con elementos de mamposterfa, 
tncluso muchos de ellos estan situados en zonas sfsmicas. Resulta por ello de fundamental 
importancia, poder evaluar su comportamiento bajo acciones estaticas y dinamicas. 

La historia del comportamiento de edificios de mamposterfa sometidos a acciones sfsmicas, 
muestra una larga lista de fallas. Sin embargo, como el mismo analisis lo demuestra, las causas de 
vulnerabihdad son evitables y las edificaciones de mamposterfa correctamente disenadas, 
detalladas y construidas han demostrado competencia sismorresistente en los sismos mas severos 
A pesar de la importancia que tiene, el analisis y la prediction de la respuesta mecanica de 
elementos estructurales de mamposterfa es todavfa un problema abierto. Se han realizado 
numerosas investigaciones para la determinacion de las relaciones constitutivas del material y 
para la defimcion de un modelo numerico efectivo, sin embargo, el problema es diffcil de abordar 
por muchas razones, entre ellas, la gran diversidad de materials y tipologfas estructurales, 
ademas, la medicion directa de las propiedades mecanicas de la mamposterfa es bastante diffcil 
por muchas razones, entre las cuales pueden mencionarse: 1) la mamposterfa mas que un material 
compuesto es una estructura y, por ello, sus propiedades son fuertemente dependientes de la 
disposicion de sus elementos y su interaccion; 2) la escala a la que deben hacerse los ensayos 
depende de la micromecamca del fenomeno, que debe ser tenido en cuenta; 3) la gran dispersion 
en los resultados experimentales, debido a la naturaleza no homogenea de los materiales. a la 
mano de obra, y a las inevitables imperfecciones geometricas, requieren definiciones especfficas 
y precisas de los procedimientos experimentales; 4) generalmente las metodologfas adoptadas por 
diferentes investigadores no son comparables; 5) las definiciones de los parametros mecanicos de 
la mamposteria y de los procedimientos experimentales para su determinacion, no son 
congruentes ya que no estan referidos a un mismo modelo micromecanico, en general, se usan 
ensayos originalmente disenados para materiales lineales elasticos. 

El objetivo de este trabajo es simular numericamente el comportamiento de elementos pianos 
de mamposterfa. Para ello, se propone una metodologfa y un modelo constitutivo que permiten 
simular el comportamiento de estos elementos, a partir de las propiedades mecanicas de los 
materiales constituyentes y de su disposicion geometrica. 

En primer lugar se presenta una breve introduccion al comportamiento de la mamposterfa, un 
resumen del estado del arte y una clasificacion de los metodos para el analisis de la mamposterfa. 

Postenormente se presenta el modelo propuesto, se indican las etapas de calculo, se describe 
el modelo numerico para materiales anisotropos utilizado, y se indica la metodologfa para la 
obtencion de los parametros mecanicos de la mamposterfa. 

Por ultimo, se muestran algunos ejemplos de aplicacion en muretes de ladrillos huecos en 
muros de mamposterfa simple con carga en el piano, en porticos rellenos con mamposterfa 

sometidos a carga lateral y en muros sometidos a cargas cfclicas, que prueban la validez del 
modelo propuesto. 



Estudio (coric'o-num6rico del comportamiento de la marnpasterfa 



3 



2 COM PORTAMI ENTO DE LA MAMPOS TERIA 

Las propiedades de la mampostena dependen de las propiedades de los ladrillos. de las 
propiedades del mortero y de la disposicion de los ladrillos y las juntas (topologfa estructural), en 
consecuencia. se trata no de un material compuesto, sino de una estructura heterogenea v 
anisotropa. 

El comportamiento de la mampostena sujeta a estados complejos de tension, esta nu\ 
influenciado por la orientacion de las juntas de mortero respecto de las cargas. La mampostena 
presenta un comportamiento marcadamente direccional, justamente debido a la influencia de las 
juntas de mortero que actuan como pianos de debilidad (Page 1981). 

Las propiedades elasticas de la mampostena, varian con la inclinacion de la junta de mortero 
respecto de la carga aplicada. Sin embargo, dicha variacion, es bastante moderada en relacion con 
la diferencia que existe en los modos de falla y en la resistencia cuando varfa la inclinacion de la 
junta. La falla no puede ser definida a traves de un criterio en terminos de tensiones principals 
en un punto. Se debe tener en cuenta tambien la orientacion de las tensiones principals respecto 
de las juntas de mortero. Dependiendo de dicha orientacion, la falla puede ocurrir solo en la junta 
o en un mecanismo combinado de rotura de juntas y ladrillos (Page 1981). 

Por otro lado, el tipo de direccionalidad en el comportamiento, depende tambien de la 
disposicion geometrica de los mampuestos en el elemento. Si se considera un bloque de 
mampostena como el de la Fig. La, se tiene un caso tfpico de isotropfa polar, en el que el eje de 
simetrfa material es perpendicular al piano del mortero. Distinto es el caso de un panel de 
mampostena, que debe ser considerado ortotropo con tres ejes de simetrfa material como se 
indica en la Fig. Lb. 




Figura 1 a) Bloque con isotropfa polar, b) Panel con ortotropfa material. 



Ensayos de compresion cfclica en paneles de mampostena (Naraine 1991, Bernardini 1978 
1979) han dado como resultado curvas similares a las obtenidas para el hormigon. En la Fig.2 se 
ha representado dicha curva para compresion normal a la junta de mortero. En la misma se pone 
en evidencia la presencia de deformaciones permanentes y degradacion de rigidez bajo carga. 

La mampostena transmite fuerzas .de compresion en forma muy efectiva. Su capacidad de 
absorcion de fuerzas de compresion esta gobernada, entre otras cosas, por la resistencia a traccion 
de los ladrillos (Page 1981), los que fallan debido a la tension transversal provocada por la 
diferente expansion lateral que tienen estos ultimos (mas rfgidos), respecto al mortero (mas 
flexible) (Page 1981). La resistencia a traccion de la mampostena es extremadamente baja. 
debido a la baja adherencia entre el mortero y los ladrillos. Para traccion normal a la junta se 
obtiene la minima resistencia (Page 1981). 



4 



P. Martin. B. Luccioni y R. Danesi 



En un panel de mamposterfa sujeto a cargas en el piano, se pueden distinguir dos modos de 
falla dependiendo de la orientacion de las tensiones principales. Ejemplos tfpicos de estos modos 
de falla se muestran en la Fig. 3. 

<j[NVmi>] 

6t 




Figura 2 Ensayo de compresion ciclica uniaxial normal a la junta. 



Para compresion uniaxial la falla ocurre en pianos normales al del panel. Dependiendo de la 
orientacion de las juntas de mortero respecto de la carga aplicada, la falla ocurre por fisuracion y 
deslizamiento en la cabeza o pie de los paneles o de una forma combinada de fisuracion de juntas 
y ladrillos (Page 1981). Para compresion uniaxial paralela al piano del mortero, la falla ocurre 
inicialmente por expansion lateral de las juntas de mortero, Fig.3.a (Page 1981). Las columnas de 
ladrillo resultantes son capaces de seguir soportando carga, pero la primera falla es tomada como 
carga de colapso y corresponde al momento en que las deformaciones laterales del panel 
aumentan bruscamente (Page 1981). 








CARGA 




Figura 3 Falla de paneles con carga en el piano, para distintas orientaciones respecto a los ejes de 

ortotropfa. 



En el caso de compresion biaxial, esos modos de falla son impedidos por la presencia de una 
segunda tension principal de compresion G r Para la mayoria de las relaciones g^g,, la falla 



Estudio te6rico-num€rico del comportamiento de la mamposterfa 



5 



ocurre en pianos paralelos al piano del panel, independientemente de la orientacion de la junta de 
mortero (Page 1981) (ver Fig.4). 

La transicion entre los dos modos de falla de las Figs. 3 y 4 se produce para valores altos de la 
relacion gJg 2 y depende de la orientacion de la junta. 

Las envolventes de falla promedio, obtenidas para cada orientacion de las juntas han sido 
representadas en la Fig.5. 

Resulta importante hacer notar que la relacion entre resistencia a compresion y traccion, no es 
la misma en las dos direcciones principales de ortotropfa. Esto es: 



( T, CJ-i 

I comp 3 comp 

G , (7, 

1 trace 3 trace 



( 1 ) 



3 MODELOS PARA MAMPOSTERIA 



La unica manera de controlar los estados lfmite de servicio, entender a fondo los mecanismos 
de falla y poder evaluar de forma confiable la seguridad estructural, es a traves de modelos 
constitutivos robustos, complementados con tecnicas numericas eficientes para la solucion de los 
sistemas de ecuaciones que resultan de la discretizacion de los problemas (Louren^o 1 996) 




Figura 4 Falla dc panel bajo carga biaxial de compresion. 



En el analisis mediante elementos finitos de estructuras de mamposterfa sin armar, el efecto de 
las juntas de mortero como causa fundamental de la anisotropfa y la no-linealidad en el 
comportamiento, ha sido tenido en cuenta con diferentes grados de refinamiento. En general, la 
representacion numerica del comportamiento se puede centrar en una micro-modelacion de cada 
uno de los componentes individuates, esto es mampuestos y mortero o en la macro-modelacion de 
la mamposterfa como un compuesto, (Rots 1991). 
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Figura 5 Ensayo de compresion biaxial para distintos dngulos (Page, 1978). 



El objetivo principal de los micro-modelos es simular adecuadamente el comportamiento de la 
mampostena a partir de las propiedades mecanicas de los componentes y de las interfaces. 
Representan un metodo adecuado para elementos estructurales pequenos en los que se quiere 
estudiar especialmente estados de tension y deformacion fuertemente heterogeneos. Los macro- 
modelos son aplicables cuando la estructura es lo suficientemente grande como para tener 
tensiones aproximadamente uniformes en macro longitudes. Claramente este tipo de enfoque esta 
mas orientado a la practica, debido a la menor demanda de tiempo y memoria, y a la facilidad 
para generar la malla de elementos finitos. 

3.1 Micro-modelos 

Dentro de los micro-modelos se pueden distinguir dos tipos de estrategias: Los micro-modelos 
detallados en los que los que los mampuestos y las juntas son representados por elementos 
continuos, mientras que la interfaz entre ambos se representa mediante elementos discontinues 
(Page 1978, Lofti 1994), y los micro-modelos simplificados, en los que los mampuestos se 
representan mediante elementos continuos, mientras que el comportamiento de las juntas y de las 
interfaces junta-mampuesto es considerado en forma conjunta, mediante elementos de interfaces 
(Rots 1991, Louren^o 1996). 

En el primer enfoque se tienen en cuenta todas las propiedades mecanicas del mortero y las 
interfaces representan pianos de deslizamiento potenciales. En el segundo enfoque cada junta, 
que consiste en mortero y dos interfaces, esta concentrada en una interfaz promedio. La 
mamposterfa es considerada entonces, como un conjunto de bloques elasticos unidos por 
superficies de fractura o deslizamiento potencial. Como no se considera el efecto Poisson del 
mortero, se pierde precision. 
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3.2 Macro-modelos 

En este tipo de modelos los mampuestos, el mortero y las interfaces son considerados de 
manera promedio en el continuo y se analiza a la mamposteria como un compuesto anisotropo y 
se establece una relacion entre tensiones promedio y deformaciones promedio. Un macro-modelo 
completo debe reproducir un material ortotropo con diferentes resistencias a traccion y a 
compresion y propiedades inelasticas en cada una de las direcciones de ortotropia. Se ha 
propuesto un numero reducido de macro-modelos especfficos para mamposteria (Anthoine 1992, 
Lourengo 1996, Martin 1997) y muchos de el los no han sido implementados debido a la 
dificultad para formular algoritmos numericos robustos que permitan tratar la ortotropia y las 
propiedades inelasticas eficientemente. 

Dentro de los macro-modelos propuestos para mamposteria, existen fundamentalmente dos 
tipos de enfoques: un primer enfoque claramente fenomenologico, en el que los parametros del 
material deben ser obtenidos de ensayos de elementos de mamposteria suficientemente grandes 
bajo estados de tension homogeneos (Dhanasekar 1985; Cuomo 1992; Contro y Nova 1982) y un 
segundo tipo de enfoque, comunmente denominado micro-macro modelo, en los que se utiliza 
algun tipo de tecnica de homogeneizacion para obtener el macro-modelo a partir de los micro- 
modelos de los componentes y su disposicion (Pietruszczak y Niu 1992; Anthoine 1995; Sulem y 
Mulhaus 1997). En muchos de estos modelos, la macro-ley constitutiva del compuesto no se 
implementa, ni siquiera se conoce, solo puede obtenerse el comportamiento del compuesto para 
un dado camino de carga. En otros casos los micro-modelos se utilizan para obtener los datos del 
material promedio para el macro-modelo (Louren^o 1996, Martin 1997). El modelo propuesto en 
este trabajo corresponde a este ultimo tipo de modelos. 



4 ESQUEMA DEL PROCESO DE CALCULO 

La idea del metodo propuesto, es simular el comportamiento de la mamposteria de ladrillos a 
partir de las propiedades mecanicas de los materiales intervinientes y de su disposicion 
geometrica, sin necesidad de discretizar cada uno de los elementos por separado. Para ello, se 
procede en dos etapas. En una primera etapa, se analiza el comportamiento de una celda basica, 
que se discretiza en forma fina. De los resultados de dicho analisis, se obtienen las propiedades 
mecanicas del conjunto. En una segunda etapa, se analiza ya el comportamiento de la estructura 
de interes, considerando a la mamposteria como un material homogeneo anisotropo cuyas 
propiedades mecanicas son las obtenidas en el preproceso. 

Tanto en el preproceso de la celda basica, como en la simulacion de la mamposteria como un 
material homogeneo con propiedades promedio, se utiliza un modelo constitutivo de plasticidad y 
dano acoplados (Luccioni 1996, Luccioni 1995). Dicho modelo, es adecuado para la simulacion 
del comportamiento de geomateriales, y permite tratar problemas de materiales compuestos 
anisotropos. 



5 MODELO CONSTITUTIVO 

Muchos materiales exhiben un comportamiento marcadamente direccional. La simulacion del 
comportamiento de este tipo de materiales es un problema de gran complejidad. Existen 
basicamente dos tipos de enfoques para simular el comportamiento elastoplastico anisotropo: 
formular funciones de fluencia para solidos anisotropos (Hill 1967) o simular el comportamiento 
del solido anisotropo real mediante un solido ficticio isotropo (Betten 1981, Betten 1988). 

La idea de trabajar con un solido isotropo, fue introducida por Betten, quien uso este tipo de 
enfoque para simular plasticidad y creep de materiales anisotropos (Betten 1981, Betten 1988). 
En trabajos anteriores (Luccioni 1993, Oiler 1993a, Oiler 1993b, Luccioni 1995) se utilizo esta 
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idea para simular el comportamiento de materiales elastoplasticos inicialmente anisotropos. En 
los mismos, se supone la existencia de dos espacios: un espacio anisotropo real y un espacio 
isotropo ficticio. El problema se resuelve siempre en el espacio isotropo ficticio, utilizando 
funciones de fluencia originalmente desarrolladas para materiales isotropos. Los tensores de 
tension en ambos espacios estan relacionados mediante una transformacion lineal. El tensor de 
transformacion contiene toda la informacion sobre la anisotropia inicial del material, es 
independiente del tipo de carga y se mantiene constante a medida que evoluciona el proceso 
elastoplastico. 

La mamposterfa de ladrillos, presenta un comportamiento bastante mas complejo que el que 
puede describirse mediante estos modelos. En este tipo de materiales, el grado de anisotropia o la 
relacion entre resistencias en distintas direcciones, depende del tipo de estado tensional: 
compresion, traccion o cualquier combinacion de tensiones y a la vez, evoluciona a lo largo del 
proceso elastoplastico. 

En el presente trabajo se utiliza una generalization (Martin 1997, Luccioni 1997) de los 
modelos antes desarrollados (Luccioni 1993, Oiler 1993a, Oiler 1993b, Luccioni 1995) que 
permite tener en cuenta las caracterfsticas de anisotropia descriptas. A1 igual que en los modelos 
anteriores, se resuelve el problema en un espacio isotropo ficticio. La diferencia, es que los 




Figura 6 Transformacion de espacios. 



tensores de tension en el espacio isotropo ficticio y en el espacio anisotropo real, no estan 
relacionados mediante una transformacion lineal, sino que el tensor de transformacion es funcion 
del tipo de estado tensional, y evoluciona con el proceso elastoplastico. 

5.1 Mapeo de espacios 

El modelo parte de la hipotesis de que existen dos espacios (Betten 1981, Betten 1988, 
Luccioni 1993, Oiler 1993a, Oiler 1993b, Luccioni 1995): a) Un espacio anisotropo real y b) Un 
espacio ficticio isotropo, ver Fig.6. El problema se resuelve en el espacio isotropo ficticio Io que 
permite utilizar modelos elastoplasticos desarrollados para materiales isotropos. 

Los tensores de tension en ambos espacios estan relacionados mediante una transformacion 
del tipo: 

T IJ = A ljkl (<J,K p )CT k] (2) 

Donde: T tJ y <J k] son los tensores de tension en los espacios (a) y (b) respectivamente y 

es un tensor de transformacion de cuarto orden, que contiene la informacion sobre la anisotropia 
de resistencia. En el caso mas general, este tensor es funcion del tipo de estado tensional y de la 
evolucion del proceso elastoplastico. En este trabajo se propone para el mismo la siguiente 
expresion: 
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4 »n=4jkl(K P ) S klnm +A‘ kt (K p ) ( I klm „-s klnw ) (3) 

Donde I klmn es el tensor identidad de orden 4. Los tensores A! tjkI (K p )y A- jkj (K p ) son los 
tensores de transformation correspondientes a procesos de traccion y compresion simple, y son 

funciones de la variable de dano plastico K p (Luccioni 1993, Oiler 1988a, Oiler 1988b), que 
hace las veces de variable de endurecimiento plastico isotropo. Para un dado estado 

termodinamico identificado a traves de la variable de dano plastico K p , los tensores 
A^ t {K p ) y se calculan a partir de los valores de las funciones tensoriales de 

endurecimiento en procesos de traccion y compresion simple: 



4 /(^) = fl (*') fL(>< p r' < w ( 0 = /; (^)A/(^r l (4) 

Donde y Jl/Son las funciones tensoriales de endurecimiento en traccion en los espacios 

isotropo ficticio y anisotropo real respectivamente, y J* y f kl son las funciones tensoriales de 

endurecimiento en compresion en los espacios isotropo ficticio y anisotropo real respectivamente. 
Todos estos tensores, dependen de la evolucion del proceso plastico, a traves de la variable de 

dano plastico K p , y deben estar referidos al sistema de coordenadas global. En el caso de 
materiales ortotropos, si se toma como referencia las direcciones principales de anisotropfa, estos 
tensores son diagonales. Los elementos de la diagonal son los valores de la funcion de 
endurecimiento en dichas direcciones. Para poder calcular A‘ Jk/ (fC p ) y A‘ k/ (fC r ) mediante las 
ecuaciones (4), se debe rotar dichos tensores diagonales al sistema de referencia global. 

El tensor S kimn tiene en cuenta el tipo de estado tensional y se calcula como: 

S, ,_=<*, A <» 

2(7 r, 



5.2 Caracteri'stica generates del modelo isotropo 

El modelo utilizado es termodinamicamente consistente y esta basado en la hipotesis de 
elasticidad desacoplada (Lubliner 1972). De acuerdo a esta hipotesis, la energia libre puede 
suponerse formada por dos partes independientes: una parte elastica y una parte plastica. 

(6) 

La ecuacion constitutiva secante, afectada por la degradacion de rigidez, se deriva de la 
desigualdad de Clasius-Duhem (Malvern 1969) y puede escribirse en terminos de la energia libre 
como sigue 



T u =m„ 



dV* 

del 



=q; l e' tl =(\-d)c° a ' kr 



(7) 



Esta ultima expresion (7), corresponde al caso mas simple de degradacion de rigidez, que 
normalmenie se conoce como dano escalar. d es la denominada variable de degradacion de 
rigidez que vana entre 0 para el material virgen y 1 para el material totalmente danado. El tensor 
constitutive secante degradado resulta: 



q; l ={\-d)q; l 



(8) 



5.3 Proceso plastico 

El umbral de comienzo de comportamiento plastico, se define mediante una funcion de 
fluencia (Luccioni 1993, Oiler 1993a, Oiler 1993b, Luccioni 1995): 



F{(j IJ \a k ) = FiT^CCk) = 0 



(9) 
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Donde: F y F designan las funciones de fluencia en los espacios anisotropo real e isotropo 
ficticio, y CC y OC son variables intemas correspondientes a dichos espacios. 

La transformacion de espacios definida por la ecuacion (2), permite utilizar funciones de 
fluencia F definidas para materiales isotropos en el espacio isotropo ficticio. Debe notarse, que 
dicho espacio, es isotropo en cuanto a umbrales de comportamiento plastico, pero no 
necesariamente respecto a otras propiedades, como la rigidez elastica, por ejemplo. 

La deformacion plastica en el espacio real, se calcula mediante la regia del flujo: 




Figura 7 Superticie de fallaen el espacio real anisotropo y en el espacio ficticio isotropo (Mohr Coulomb). 

Donde G es la funcion de potencial definida en el espacio de tensiones reales. Bn lugar de 
trabajar con esta funcion, que debe ser anisotropa, se puede trabajar con una funcion de potencial 

G , definida en el espacio isotropo ficticio tal que: 



G(a 0 ,a k ) = G( x, r a k ) 

La ecuacion (10) se escribe entonces: 



e p 

c 'j 



= Xf- = X-f-^ L = X^-H khJ =XR IJ 
da v dr k , oOy ox k , 



Donde: 



tt ki _ ^ (A klnw^ mn ) — A • 



dA k lmn ~ 

d °mn 



De manera analoga, la ley de evolucion de las variables intemas se escribe como. 

-~z i I.m dG ^ J dG dt kl _ i u m v 



(li) 



( 12 ) 



(13) 



(14) 



5.4 Proceso de dafio 

El lfmite entre el comportamiento no danado, y el comportamiento daiiado, se describe a 
traves de una funcion de dafio, similar a la funcion de fluencia en plasticidad. 
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G D {<y ip p k ) = G D {r ij ,p k ) = T{T iJ )-f e {P k )< 0 ( 15 ) 

Donde: G D y G D designan las funciones de dano en los espacios anisotropo real e isotropo 
ficticio, y J3 k y /3 k son variables internas de dano correspondientes a dichos espacios, t(T t/ )c s 

una funcion escalar de argumento tensorial denominada tension equivalente y / c .(/? A )es el 
umbral de dano equivalente. 




Figura 8 Superficie de falla en el espacio real anisot. y en el espacio ficticio isotropo (Von Mises). 




Figura 9 Superficie limite para distintas orientaciones de los ejes de ortotropfa (Mohr Coulomb). 
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5.5 Proceso de dano y plasticidad acoplados 



La respuesta elasto-plastica y de dano acoplada, se obtiene de la solucion simultanea del 
siguiente sistema de ecuaciones de consistencia en el espacio isotropo ficticio (Luccioni 1993, 
Luccioni 1996): 



F(T,j,a k ) = 0 

'G D \x ip P k ) = 0 



(16) 



5.6 Ejemplos ilustrativos 

En las Figs. 7 y 8 se ilustra como, partiendo del criterio de Mohr Coulomb o de Von Mises, 
respectivamente, en el espacio ficticio isotropo, y utilizando la transformation de espacios 
descripta, se puede lograr una funcion de fluencia anisotropa en el espacio real, donde la relacion 
de resistencias es distinta en ambas de direcciones. En el caso del ejemplo desarrollado, la 
relacion entre resistencias a compresion en la direccion 1 es 2.4 y para la direccion 3 es 4.8. 

Si los ejes de ortotropfa del material, estan girados respecto de los ejes globales, un cierto 
angulo, es necesario rotar los tensores de transformacion, esto es, expresarlos en tuncion de las 
coordenadas globales. En este caso, el tensor de transformacion ya deja de ser diagonal y pasa a 
ser un tensor completo. A modo de ejemplo, se muestra en la Fig.9, la superlicie If mite en el 
espacio de tensiones reales, para distintas inclinaciones de los ejes de ortotropfa respecto a los 
ejes globales de referenda. 

5.7 Datos del modelo 

En el caso de un material ortotropo general, como el descripto, el numero de datos necesarios 
para reproducir su comportamiento mecanico, en cualquier direccion, es bastante mayor que en el 
caso de un material isotropo. 

La mayorfa de los datos requeridos para el modelo que se presenta, pueden ser obtenidos a 
partir de las curvas de respuesta tension - deformacion, para traccion y compresion simple, en 
cada una de las direcciones principales de ortotropfa. 

1 . Parametros elasticos. (Para todas las direcciones de ortotropfa) 

1.1. Modulos de Young. 

1 .2. Modulos de elasticidad transversal o Modulos de Poisson. 

2. Valores de resistencia iniciales a traccion y compresion en cada una de las direcciones de 
ortotropfa. 

3. Curvas de endurecimiento en las direcciones de ortotropfa del material. 

Los ensayos que deben realizarse para determinar los parametros mecanicos de la 
mamposterfa son complejos y delicados, y esto los hace muy costosos. Por ello, resulta atractivo, 
poder obtener los parametros basicos del conjunto, a partir de los parametros mecanicos de los 
materiales componentes y de su disposition. 



6 PREPROCESO DE UNA CELDA BASICA 



Con el objeto de obtener los datos para el macro-modelo descripto, se afsla una portion de la 
mamposterfa, denominada celda basica (Fig. 10), que es representativa de la mamposterfa 
utilizada, y se la analiza en detalle en un esquema de elementos finitos, discretizando en forma 
fina cada material por separado. 
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A 





Figura 10 Panel con distintas configuraciones de celda basica y discretizacidn dc la celda basica. 



Para determinar los parametros en forma adecuada, es necesario como minimo, el estudio del 
comportamiento de la celda bajo cuatro estados de carga: 

1. - Traccion en la direccion paralela a la junta horizontal 

2. - Compresion en la direccion paralela a la junta horizontal 

3. “ Traccion en la direccion perpendicular a la junta horizontal 

4. - Compresion en la direccion perpendicular a la junta horizontal 

El metodo es aplicable, solo si el tarnaho de la celda es pequeno respecto de la estructura, y las 
propiedades geometricas y mecanicas del panel se repiten periodicamente. 

6.1 Parametros mecanicos deducidos 

Del analisis de la celda basica se pueden obtener los siguientes datos para el macro-modelo: 

1 . Modulos de Elasticidad y Modulos de Poisson en las direcciones principales de anisotropfa. 

2. Lfmites elasticos del conjunto para traccion y compresion uniaxial en cada una de las 
direcciones principales de anisotropfa. 

3. Funciones de endurecimiento de la celda para cada direccion en traccion y compresion, en 
terminos de la variable de dano plastico. 

En la Fig.l 1 se muestran los resultados obtenidos del analisis de una celda basica. En ellos se 
puede observar el comportamiento diferenciado en distintas direcciones. 




Figura 1 1 Resultados del analisis de una celda basica. 



7 EJEMPLOS DE APLICACION 

El esquema de calculo y el modelo descripto, fueron implementados en un programa de 
elementos finitos, con el cual se resolvieron los ejemplos que se presentan a continuacion. 
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7.1 Analisis de muretes de ladrillos ceramicos huecos 



En el Laboratorio de Estructuras de la Universidad 
de una empresa del medio, un estudio sobre 10 pilas a 



Nacional de Tucuman, se ejecuto, a pedido 
compresion uniaxial y sobre 10 muretes a 




Figura 12 Esqucma de la instrumcntacidn de la probeta. 





Probeta o ensayor 
Losa Reactive 



Portico de Cargo 



Gato Hidraulico 



Celdo de 
Cargo 



Figura 13 Dispositivo de ensayo. 



compresion diagonal (Martin 1995a, Martin 1995b). Los mampuestos utilizados eran ladr.llos 
ceramicos huecos de 13 cm de espesor, obtenidos por extrusion. Se instrument la ultima probeta 
de cada serie de ensayos, a fin de determinar el estado tensional real, y contrastar estos resultados 
con los obtenidos mediante los modelos teorico - numericos desarrollados. 
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Figura 14 Comparac<on dc resultados experiment ales y numcricos. 

Para la instrumen acion, se utilizaron extensometros electricos adheridos a la superfic e ■ 
mortero y del ladrillo en distint is posiciones y direcciones como se indica en la Fig I . La 
instrumentacion, se completo con la colocacion de una celda de carga de 1000 kN ae < tp.cidad 
maxima, y con un sensot de desplazamiento el que fue contrastado con un flexfme a durante ,-l 
ensayo. 

En la rig. 1 3 se ^uestra e! esquema del dispositivo de carga utilizado en los en o de 
compresnn y impresion diagonal re' pectivamente. 

Estos ens; yos >e modelaron numericamente, discretizando cada uno de los mater, ie 
interviniente' y la nterfaz por separrdo, co i la malla de elementos finitos indicada en la Fig. 
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En la Fig. 14, se resumen algunos de los resultados de los ensayos y su comparacion con los 
valores obtenidos numericamente. Los graficos a) a e) corresponden a curvas carga-deformacion 
en 




Figura 15 Esquema del panel y mallas de elemcntos finitos. 



algunos de los extensometros dispuestos y el grafico f) corresponde a la curva carga - 
desplazamiento del extremo superior. De la observation de esta> curvas, se puede concluir que, 
en terminos generates, y a pesar de la variabilidad de las propiedad^; de los materia* . y de la 
complejidad del ladrillo para ser considerado un elemento piano, el modelo numerico reproduce 
adecuadamente los resultados experimentales. 

La carga ultima obtenida numericamente (420 kN) es muy cercana a la experimen ' y la 
curva carga - desplazamiento del extremo superior, en la Fig. 1 4-f, muesLa tambien muy buena 
coincidencia entre resultados numericos y experimentales. 

La pieza fallo con la formacion de fisuras verticales siguiendo la lfnea de los extensometros de 
la Fig.l4-b, mientras que en el ensayo se formo una gran fisura en el centro del especimen. En la 
mayorfa de los especimenes de la serie, la fisuracion fue del tipo de la observada numericamente. 
A1 tratarse de materiales con tanta variabilidad y fragilidad, en los que cualquiei [ equena falla 
puede inducir la formacion de la fisura en un lado u otro del especimen, las diferencias 
encontradas son totalmente justificables. 

Debe notarse tambien, que la diferencia entre el modulo de Poisson de los ludrillos, y el 
correspondiente al mortero de las juntas, genera tensiones de traccion en la direccion horizontal 
en los ladrillos. Dichas tensiones son las responsables de las fisuras verticales que conducen a la 
rotura de ladrillos. 

7.2 Analisis de un panel de mampostena maciza 

En el siguiente ejemplo se contrasta el modelo numerico propuesto con los resultados 
experimentales obtenidos por Page (Page 1978), en un muro de mampostena maciza, sometido a 
un estado de carga uniaxial uniforme. 

El analisis del panel se realizo utilizando dos enfoques distintos, primero con un macro- 
modelo y luego con un nucro-modelo En un primer paso, una vez obtenidos los parametros 
mecanicos de la mampostena mediante el analisis de una celda basica, se procedio a simular el 
comportamiento mecanico del panel, bajo distintos estados de carga, considerando a la 
mampostena como un material homogeneo. El analisis se realizo utilizando una discretizacion 
gruesa, como la que se muestra en Fig. 15 (seis elementos isoparametricos ue 4 nodos con 2x2 
puntos de integracion), 

Posteriormente, el mismo panel se discretizo de una forma fina. considerando juntas de 
mortero y mampuestos como elementos independientes. Se utilizaron elementos isoparametria \s 
de A nodos con 2x2 puntos de integracion, ver Fig. 15. 

En la Fig. 17, se muestran los resultados obtenidos en ambos tipos de analisis, y su contraste 
con los resultados experimentales, para el caso del panel sometido a una carga cuasiestatica de 
cernpresion perpendicular a la junta horizontal. 
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Se realizo el mismo analisis para carga paralela a la junta horizontal. En este caso solo se 
compararon los resultados provenientes de una discretizacion fina, y los de una discretizacion 
gruesa con el material homogeneizado (Fig. 17), ya que no se dispoman de resultados 
experimentales. 

Por ultimo, se analizo el mismo panel sometido primeramente a una carga de compresion 
normal a la junta y luego a una carga horizontal creciente en su extremo superior que le producia 
solicitaciones de tipo cortante. Se considero una discretizacion fina, distinguiendo cada uno de 
los materiales, y dos discretizaciones gruesas distintas, considerando al material con sus 
propiedades promedio.' Los resultados obtenidos se presentan en la Fig. 18. En este caso, se 
observa tambien una buena concordancia, entre las distintas formas de analisis, teniendo presente 
que el estado tensional es mucho mas complejo, que los analizados previamente. 

El ejemplo descripto se resolvio con las distintas configuraciones de celdas basicas indicadas 
en la Fig. 19, observandose que, si bien las diferencias obtenidas no son grandes, a medida que el 
tamaho de la celda aumenta, la respuesta es menos rigida, y tiende a la solucion obtenida 
discretizando todo el panel en forma fina. 

7.3 Portico de hormigon relleno con mampostena 

Mehrabi et al (Mehrabi 1995), realizaron una serie de ensayos sobre porticos de hormigon 
rellenos con mampostena. Los paneles se sometieron a una pre-compresion normal a la junta 
horizontal y posteriormente se los solicito con carga lateral. En la Fig. 20 se esquematizan los 
paneles y dispositivos de carga. 

La malla de elementos finitos se indica en la Fig.21. Para los elementos de mampostena, se 
utilizo un macro-modelo, con las propiedades mecanicas obtenidas del preproceso de una celda 
basica. Para tener en cuenta la armadura del portico de hormigon armado se aplico, para el 
mismo, un modelo de compuesto basado en la teona de mezclas (Truesdell 1960; Luccioni 1993) 

En la Fig.22 se muestran las curvas carga-desplazamiento lateral del extremo superior del 
panel obtenidas numerica y experimentalmente. En la misma se observa que el modelo reproduce 
ajustadamente el comportamiento observado experimentalmente. 

Oy 

[MPa] 




Figura 16 Comparacion de resultados experimentales, y numcricos. 
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a* 




Figura 17 Comparacion de resultados numericos con ierando una discretizacidn fina, y una discretization 

gruesa. 




desplazamiento lateral 
d[m] 

Figura 18 Muro de mamposterfa con carga lateral. 



rm 






celda 1 celda 2 celda 3 celda 4 celda 5 
Figura 19 Distintas configuraciones de celdas estudiadas. 



7.4 Panel de mamposteria bajo carga dclita 

En este ejemplo, se presenta la aplicacion del macro-modelo descripto, a la simulacion del 
comportamiento de paneles de mamposteria ensayados por Tomazevic et al (Tomazevic 1997) 
bajo carga cfclica. Se trata de modelos a escala reducida d^ paneles de mamposteria de bloques 
con confinamiento, cuyas dimensiones pueden verse en la Fig.27. Los paneles tueron sometidos a 
una carga de compresion vertical y luego a desplazamientos horizontales cfclicos. 




Estudio te6rico-num6nco del comportamicnto de la mamposteria 



19 



Dispositivos de c:>licaci6n de cargos 




Figura 20 Esquema dei especimen y del dispositivo de ensayo. 




Figura 21 Mallu de elememos 1 utos utilizada. 




Figura 22 Comparacion de resultados ohtenidos numc; amente y experimentalmente. 
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Figura 23 Esquema y dimensiones del modelo de mampostena de bloques a escala reducido. 




Figura 24 Curvas carga-desplazamiento horizontal del extremo superior. 

En las Fig.24 se han representado la curva carga-desplazamiento horizontal del extremo 
superior obtenida numericamente, su envolvente y la envolvente de desplazamiento medida 
experimental mente. De las mismas se deduce que el modelo reproduce adecuadamente el 
comportamiento de la mampostena bajo cargas de tipo cfclico. 



8 CONCLUSIONES 

Se ha presenlado un modelo que permite reproducir el comportamiento de una amplia gama de 
materiales ortotropos, entre ellos la mampostena de ladrillos. El mismo, considera a la 
anisotropia, como un fenomeno evolutivo en funcion del proceso elastoplastico y del estado 
tensional del punto. La condition de fluencia es verificada en el espacio ficticio isotropo. por lo 
que no es necesaria la definicion de criterios de falla y de potencial en el espacio anisotropo. La 
transformation de espacios no lineal, que es necesaria para la definicion del modelo solo precisa 
de algunas transformaciones y derivadas adicionales respecto del caso isotropo. 
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Se ha implementado una metodologia sencilla y eficiente para la determinacion de los 
parametros mecanicos a ser introducidos en un modelo con no-linealidad material y anisotropia, 
tendiente a simular el comportamiento de paneles de mampostena. Esta metodologia permite 
reproducir el comportamiento de la mampostena a partir de las propiedades mecanicas de los 
materiales intervinientes y de su disposicion, sin necesidad de recurrir a ensayos experimentales. 

Del analisis de los ejemplos numericos desarrollados, se concluye que el modelo reproduce 
adecuadamente el comportamiento de la mampostena, para solicitaciones coincidentes con las 
direcciones principals de anisotropfa. Para solicitaciones combinadas es preciso la definition de 
un criterio de falla mas ajustado. 
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RESUMEN 

Santa Fe de Bogota, capital de Colombia, es una urbe de 7 millones de habitantes que se 
encuentra localizada en una zona de moderada amenaza sfsmica de Sur America. La ciudad tienc 19 
alcaldfas menores con notables diferencias ffsicas, economicas y sociales. En conjunto, el riesgo 
sfsmico de Santa Fe de Bogota es alto, pero las estimaciones comparativas de riesgo de sus 
localidades son muy distintas debido a sus diferencias sociales, economicas y de resiliencia. Aunque 
la ciudad cuenta con una microzonificacion sfsmica detallada y con escenarios de perdidas por 
terremoto estimados, ha sido necesario analizar otras variables importantes de cada alcaldfa menor 
para estimar el riesgo sfsmico real en cada una. Se desarrollo un estudio comparativo y holfstico 
para incluir aspectos sociales y economicos adicionales a las variables sismologicas y de ingenierfa 
obtenidas del estudio detallado de microzonificacion. Este artfculo presenta el enloque 
multidisciplinario usado, la rnetodologfa. los resultados y como la administracion de la ciudad 
podrfa usarlos, desde la perspectiva de la planificacion sectorial, para promover medidas de 
prevencion de acuerdo con la estimacion del riesgo de cada alcaldfa menor. 



ABSTRACT 

Bogota, the Capitol City of Colombia, is a 7 million people urban centre placed in a moderate 
seismic hazard prone area of South America. The city has 19 districts with significant physical, 
economical, and social differences. As a whole, the seismic risk of Bogota is high, but the 
comparative risk results of its districts are very different depending on social, economic, and 
resilience differences. Although the city has evaluated detailed seismic microzonation and 
earthquake loss scenarios, it was necessary to analyse other important variables at district level to 
estimate the real seismic risk of each one. A comparative and holistic study was developed to 
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include social, economical issues besides the seismological and engineering variables obtained o! 
the detailed microzoning study. This paper presents the multidisciplinary approach used, the 
methodology, the results, and how the city administration might use them, from the sectoral 
planning perspective, to promote mitigation measures according to the risk estimation of each 
district. 

1 ENFOQUE HOLISTICO 

En los ultimos anos, desde la perspectiva de los desastres naturales, el riesgo se ha intentado 
dimensionar, para efectos de la gestion, como las posibles consecuencias economicas, sociales y 
ambientales que pueden ocurrir en un lugar en un tiempo determinado. Sin embargo, la 
conceptualizacion del riesgo no ha sido integral sino fragmentada, de acuerdo con el enfoque de 
cada disciplina involucrada en su valoracion. Para estimar el riesgo de acuerdo con su definicion 
es necesario tener en cuenta, desde un punto de vista multidisciplinario, no solamente el dano 
ffsico esperado, las victimas o perdidas economicas equivalentes, sino tambien factores sociales. 
organizacionales e institucionales, relacionados con el desarrollo de las comunidades. En la 
escala urbana, por ejemplo, la vulnerabilidad como factor intemo de riesgo, debe relacionarse no 
solamente con la exposicion del contexto material o su susceptibilidad fisica de ser afectado, sine 
tambien con las fragilidades sociales y la falta de resiliencia de la comunidad expuesta. La falta 
de organizacion institucional y comunitaria, las debilidades en los preparativos para la atencion 
de emergencias, la inestabilidad politica y la falta de salud economica de un area geografica 
contribuyen a tener un mayor riesgo. Por lo tanto, las consecuencias potenciales no solo esian 
relacionadas con el impacto del evento, sino tambien con la capacidad para soportar el impacto y 
las implicaciones del impacto en el area geografica considerada. 

La evaluacion del riesgo puede llevarse a cabo mediante la siguiente formulacion general: Una 
vez conocida la amenaza o peligrosidad H h entendida como la probabilidad que se presente un 
evento con una intensidad mayor o igual a # durante un periodo de exposicion /, y conocida la 
vulnerabilidad V e , entendida como la predisposicion intrinseca de los elementos expuestos a ser 
afectados o de ser susceptibles a sufrir una perdida ante la ocurrencia de un evento con una 
intensidad /, el riesgo R ie puede entenderse como la probabilidad de que se presente una perdida 
sobre el elemento e, como consecuencia de la ocurrencia de un evento con una intensidad mayor 
o igual a /, 



Rie\t=(Hi. V e ) 1 / ( 1 ) 

es decir, la probabilidad de exceder unas consecuencias sociales y economicas durante un periodo 
de tiempo t dado (Cardona, 1 986). 

Aquf se propone una conceptualizacion del riesgo desde la perspectiva de considerar no 
solamente variables sfsmicas y estructurales, sino tambien variables economicas, sociales, de 
capacidad de respuesta o recuperacion posterremoto, o resiliencia, que permitan orientar de 
manera efectiva las decisiones de mitigacion del riesgo. Una estimacion de este tipo podria 
considerarse holistica, es decir integral o completa del riesgo. Para el efecto, se parte de las 
estimaciones de perdidas o escenarios urbanos de danos por sismo, debido a que son el resultado 
de la convolucion de la amenaza sismica, o microzonificacion de la ciudad, y la vulnerabilidad 
ffsica de las edificaciones y de la infraestructura; aspectos a partir de los cuales se define un 
indice de riesgo ffsico o “hard”. Posteriormente, se valora un indice de riesgo del contexto o 
-soft”, resultado de la estimacion del descriptor de amenaza sismica relativa y su convolucion 
con el descriptor de vulnerabilidad del contexto, que se basa en indicadores de exposicion, 
fragilidad social y resiliencia relativa de las unidades de analisis que conforman el centro urbano. 
Dichas unidades de analisis pueden ser alcaldias menores, distritos o localidades con alguna 
autonomfa administrativa, que faciliten a la autoridad administrativa de la ciudad a llevar a cabo 
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la gestion del riesgo utilizando los resultados de la estimacion holfstica y multidisciplinaria del 
riesgo sfsmico urbano. 



2 METODOLOGIA 

El procedimienfo propuesto para estimacion holfstica y relativa del riesgo sfsmico urbano 
parte de la dentifieacion de unas unidades de analisis, A, que son las areas para las cuales se 
determina el fndice de riesgo sfsmico total 5 Rt a , el cual esta expresado por* 

IRt a = IRh* . 8irh a + IRs a . Sirs* (2) 

donde IRh a , es el fndice de riesgo sfsmico duro o ffsico {hard), el cual esta basado en 
descriptores ob enidos de los estimaiivos de las perdidas potenciales urbanas que pueden causar 
futuros terremotos; IRs*, es el un ndice de riesgo sfsmico del contexto o blando {soft), obtenido 
del produefo escalade de los descriptores de amenaza sfsmica y de vulnerabiiidad del contexto, y 
8irh a , 5 irs a son los factores de participacion de cada indice para cada area de analisis A; por su 
parte 



JRh* = Z/XiR/. 8iR| ( 3 ) 

donde Xir es el valor de cada indicador / obtenido de la informacion de los escenarios de 
perdidas y 8'R,- el factor de participacion de cada indicador i, para cada are, de analisis A; y 

IRs a = a{{Hs k - /3){Vs k - /3) + 0) (4) 

siendo Hs A . el descriptor de amenaza sfsmica del contexto, Vs*, el descriptor de vulnerabiiidad del 
contexto, y ay [5 constantes de visualization relativas a la media y la desviacion estandar de los 
valores, que se mencionan adelante en la tecnica de escalamiento. A su vez 

Hs* = Z/Xh/.Sh, (5) 

tendo Xh/ e! valor de los indicadores i obtenido del estudio de microzonificacion sfsmica urbana 
y 5 h/ el facto : de participacion de cada indicador i, para cada area de analisis A: 

Vs* = Eva.8^ a +Fva.8f a +Rv a .8r a (6) 

donde, Ev a Fv A . Rv a son indicadores de exposicion, fragilidad social y falta de resiliencia, y Se a , 
8f a 8r a son sus factores de participacion para cada area de analisis A, lo que es equivalente a 

Vs A = (X / Xe i . 8e 8e a + (2 i Xf/ . 8f/) 8f a + (Z Xr,- . Sr,) 8r a (7) 

siendo Xe/, Xr 7 , Xh/, los valores de los indicadores i que componen la exposicion, fragilidad social y 
falta de resiliencia y 8 e/, 8f/, 8r,-, a participacion de cada indicador /, para cada area de analisis A, 
respectivamente. 

indices, descriptore 1 factores e indicadores se deben definir con base en informacion 
disponible pa^a todas las unidades de analisis. Conceptualmente deben reflcjar, de la manera mas 
directa posible, lo que se quiere valorar y debe evitarse el uso simultaneo de variables o indicadores 
que expresen aproximadamente el mismo aspecto. La figura 1 indica la composicion de los indices 
de riesgo de acuerdo con la nomenclatura mencionada. 
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X^jOpcratividad en cmergencias 

f igura 1 . Indices oe nesgo .eiativo para las Alcaldias Menores de Santa Fe dc Logota con los factores de participacion utilizados. 
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1.1 ESCALAMIENTO 

Antes de integrar mediante una combinacion lineal los indicadores en descriptores y a su vez los 
descriptores en factores e indices, estos deben ser escalados en unidades compatibles que permitan 
hacer analisis relativos conmensurables. El area el espacio publico para la atencion masiva de 
personas y el personal de rescate, por ejemplo, no pueden relacionarse en forma directa, porque en 
el primero se utilizan m 2 y en el segundo personas. La tecnica adoptada para este caso es escalar con 
respecto a la media X, y la desviacion estandar S /, asf: 

x, t -(x,-/?s,) 6 -x, t +(x,+/?s,) (8) 

,k as, ,k as, 

donde las X',* y X,* on el valor escalado y el dato crudo, para la localidad k y el indicador /, y ay fi 
son constantes de visualizacion relacionadas con la media y la desviacion estandar; asf si ores 0.8 y 
P es 4.0 los valores de los indicadores se escalan con respecto a la media (0.0) mas 5.0 (/3/a) y la 
desviacion estandar ( 1 .0) se expande 1 .25 (Mot). Estas constantes se deben mantener durante todo el 
procedimiento de escalamiento. La ecuacion cambia de signo cuando el indicador es inverso al 
factor que se valora, como es el caso de los indicadores de resiliencia, los cuales son inversos a la 
vulnerabilidad del contexto. 

Expresar el resultado como una combinacion lineal implica que no existe interaccion entre los 
indicadores o entre los indicadores y los factores de participacion utilizados para la ponderacion. 
No obstante, casi todos los indices de este tipo, desarrollados hasta la fecha, usan un enfoque 
basado en una combinacion lineal y la busqueda de otros enfoques ha permitido concluir que la 
combinacion lineal es aceptable, si se tienen en cuenta las incertidumbres e imprecisiones 
inherentes de los datos. Por otra parte, hacer de un fndice riesgo una funcion no lineal podrfa 
llegar a ser mas adecuado, pero no es claro que tipo de funcion podrfa ser y su complejidad 
asociada hace que el enfoque termine por ser inconsistente con el objetivo de buscar una 
metodologfa sencilla para obtener un fndice facil de evaluar y de facil entendimiento. 

1.2 PONDERACION 

Una vez los indicadores son conmensurables, se debe realizar su ponderacion. El factor de 
participacion de cada indicador ilustra que tan importante es el indicador con respecto a los 
demas en la determinacion de un factor componente o del mismo fndice de riesgo sfsmico. Un 
fndice de esta naturaleza debe tratar de capturar el conocimiento colectivo de todos los expertos 
para definir su valor, concebido por un grupo como un todo. La ponderacion es correcta siempre 
que cumpla con esa meta. Una variedad'de tecnicas de ponderacion ha sido usada para la 
construccion de otros fndices compuestos, sin embargo todos aquellos propuestos con base en 
tecnicas estadfsticas requieren, o que las variables dependientes puedan ser medidas directamente 
(por ejemplo, regresion), o que los indicadores esten bien correlacionados (por ejemplo, analisis 
de componentes principales). Dado que el riesgo no puede ser medido directamente y los 
indicadores no estan bien correlacionados, cualquiera de las tecnicas con esos prerrequisitos no 
es mas efectiva. La unica opcion que queda para la ponderacion es la evaluacion subjetiva por 
parte de expertos. Aunque esta tecnica parece ser indeseable por algunos, por la falta de una base 
explfcitamente replicable, es la unica opcion factible y razonable para capturar el criterio, la 
experiencia y el juicio de expertos (Davidson, 1997). Los factores de participacion se pueden 
obtener mediante redes neuronales que lleven a cabo la defuzzification de las variables 
formuladas a traves de conjuntos difusos (fuzzy sets). Este tipo de tecnica facilita, mediante 
funciones de pertenencia, el manejo de apreciaciones cualitativas de expertos y variables 
lingufsticas de valoracion. La descripcion uso detallado de esta tecnica esta fuera del alcance de 
este artfculo. 
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2 ESTIMACION RELATIVA DEL RIESGO SISMICO 

Mediante la metodologia antes descrita, para un centre urbano conformado por un conjunto de 
areas suburbanas, tales como alcaldfas menores, distritos, localidades o comunas, la determinacibn 
del riesgo sfsmico relativo para cada una de las unidades de analisis se obtiene de estimar el indice 
de riesgo sfsmico ffsico y el fndice de riesgo sfsmico del contexto. La agregacion de estos dos 
indices permite llegar a una valoracion holfstica o integral del riesgo. 

En cada ciudad existen indicadores o evaluaciones de orden sfsmico, poblacional. estadfstico. 
etc. que pueden utilizarse para determinar de manera aspectos que representen o reflejen las 
variables que se desean valorar. Aquf se presenta la estimacion de indices, descriptores. indicadores 
y componentes para las 19 alcaldfas menores de Santa Fe de Bogota. 

2.1 INDICE DE RIESGO FISICO 

Para la determinacion del fndice de riesgo ffsico (hard), en el caso de Santa Fe de Bogota se 
definieron los siguientes descriptores, los cuales fueron normalizados de acuerdo con el area de cada 
alcaldfa local: 

1. Area destruida por sismo, XiRi : Definida como el area probable de destruccibn en nr de las 
zonas construidas, estimada segun la metodologia ATC-13/FEMA, modificada por la 
Universidad de los Andes, en el marco de un escenario sfsmico hipotetico para la ciudad. 

2. Numero defallecidos, Xir 2 ; Numero de heridos , Xir, : Definidos como el numero probable de 
muertos y heridos, utilizando estimativos similares a los propuestos por Whitman et a! (1973) y 
ajustados en el ATC-13, en el marco de un escenario sfsmico hipotetico para la ciudad. 

3 . Roturas de la red de acueducto, X1R4; Roturas de la red de gas Xir 5 ; Longilud caida de redes 
electricas X® 6 : Definidos como el numero probable de roturas que se presentarfan en las redes 
de acueducto y gas y la longitud cafda de la red electrica, estimados segun la metodologia de la 
ATC-13, en el marco de un escenario sfsmico hipotetico. 

4. Numero de centrales telefonicas afectadas, XIR 7 ; Numero de subestaciones electricas 
afectadas, XIR g : Definidos como el numero de las centrales telefonicas y subestaciones 
electricas con una alta vulnerabilidad sfsmica de acuerdo con la evaluacion simplificada 
desarrollada con base en el ATC-21 y otros parametros adicionales. 

La figura 2 ilustra los resultados del fndice de riesgo ffsico, obtenidos para las diez y nueve 
alcaldfas menores de la ciudad. 




Figura 2 Indice de Riesgo Ffsico basado en los escenarios de perdidas, 
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Los descriptores que componen el indice de riesgo ffsico son valores estimados de la 
convolucion de la amenaza sfsmica, obtenida de la microzonificacion de la ciudad y de la 
vulnerabilidad ffsica de las edificaciones y la infraestructura de servicios publicos. Una description 
resumida de dicha evaluation se encuentra en el Apendice A. 

2.2 INDICE DE RIESGO DEL CONTEXTO 

El fndice de riesgo del contexto (soft), para Santa Fe de Bogota, se definio como el resultado del 
producto escalado del factor de amenaza sfsmica y del factor de vulnerabilidad del contexto para 
cada alcaldfa local. Este producto se debe a que la amenaza y la vulnerabilidad son mutuamente 
condicionantes y concomitantes para que exista el riesgo; lo que equivale a la convolucion de los 
descriptores que lo componen A continuation se definen estos descriptores y se relacionan sus 
indicadores componentes: 

1. Amenaza sismica del contexto Hs*: Definida como una agregacion de valores que expresan el 
nivel de peligro sismico que se presenta en el area que cubre la alcaldfa menor, caracterizada 
por particularidades relativas a la action sfsmica. Para su determination, en este caso, se 
definieron los siguientes indicadores: 

1.1 Aceleracion espectral , Xhi: Definida como el promedio ponderado del valor espectral de 
aceleracion sfsmica para perfodos T entre 0.2s y 0.5s de las zonas de microzonificacion sfsmica 
que tienen influencia en el area de la localidad considerada. 

1.2 Area de suelos blandos, Xh 2 : Definida como el porcentaje de area de la alcaldfa menor que es 
susceptible a la amplification sfsmica por las caracterfsticas dinamicas del suelo. 

1.3 Area con potencial de licuacion , Xh 3 : Definida como el porcentaje de area de la alcaldfa menor 
que presenta suelos arenosos no consolidados y saturados con alto potencialde licuacion en caso 
de sfsmos moderados o fuertes. 

1.4 Area con susceptibilidad al deslizamiento, Xh 4 : Definida como el porcentaje de area de la 
alcaldfa menor que presenta zonas con potencial inestabilidad de laderas en caso de sfsmos 
moderados o fuertes. 

La figura 3 ilustra el descriptor de amenaza sfsmica del contexto para las 19 alcaldfas menores de 
la ciudad. 
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Figura 3 Descriptor de la amenaza sfsmica del contexto. 



2- Virfnmtbilulad sismica de! contexto. Vs*: Definida como la agregacion de valores que expresan 
aspectos de exposicion poblacional, economica y estrategica, ausencia de desarrollo economico 
y social, debilidades para absorber el impacto, deficiencias en la gestion institucional y falta de 
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capacidad para la respuesta en caso de emergencia, que reflejan y diferencian la vulnerabilidad 
global de cada alcaldfa menor. Para su deterrninacion, se definieron tres indicadores: 

2.1 Exposition , Ev A : Definida como el volumen normalizado de poblacion, edificaciones y bienes 
economicos y estrategicos expuestos a la accion sfsmica en el area que cubre la alcaldfa menor. 
Para su deterrninacion, en este caso, se definieron los siguientes componentes: 

2.1.1 Poblacion , Xei: Definida como el numero promedio de habitantes expuestos en miles en el 
area que cubre la jurisdiccion de la alcaldfa menor. 

2.1.2 Densidad poblacional, Xe 2 : Definida como el numero habitantes dividido por el area 
construida, lo que refleja el grado de concentration y congestion de personas en el area de 
la alcaldfa menor. 

2.1.3 Area construida , Xe 3 : Definida como el area normalizada de edificacion y desarrollo urbano 
en el area que cubre la alcaldfa menor. 

2.1.4 Area industrial , Xe 4 : Definida como el area normalizada de zonas industriales, fabricas o 
empresas en el area de la alcaldfa menor, lo cual refleja valores y dependencia economica. 

2.1.5 Area institutional Xe 5 : Definida como el area normalizada de zonas de uso institucional o 
gubernamental en la jurisdiccion de la alcaldfa menor, lo cual refleja valores estrategicos y 
dependencia polftica. 

2.2 Frapilidad social Fv*: Definida como el inverso del deficit de desarrollo, economico y social 
caracterizado por las condiciones de marginalidad, bajo nivel de salubridad, delincuencia y 
necesidades basicas insatisfechas de la poblacion ubicada en el area que cubre la alcaldfa 
menor. Para su deterrninacion, en este caso, se definieron los siguientes componentes: 

2.2.1 Area de barrios marginales , XF,: Definida como el area normalizada de asentamientos 
humanos ilegales con deficiencia de servicios publicos y de baja estratificacion 
socioeconomica en el area que cubre la jurisdiccion de la alcaldfa menor. 

2.2.2 Mortal idad, Xf 2 : Definida como la tasa o numero de personas que mueren anualmente por 
causa natural por cada mil habitantes en la jurisdiccion de la alcaldfa menor. 

2.2.3 Delincuencia , Xf 3 : Definida como la tasa o numero de delitos anuales por cada mil habitantes 

en el area de la alcaldfa menor, el cual representa deterioro social en la zona. 

2.2.4 Disparidad social Xf 4 : Definida como el nivel de necesidades basicas insatisfechas y de 
desarrollo humano relativo de la alcaldfa menor. 

2.3 Fait a de resilience Rv*: Definida como el inverso de la capacidad economica, social e 
institucional (resiliencia), que representa la debilidad para absorber el impacto de una crisis, la 
falta de capacidad para responder en caso de emergencia y las deficiencias en la gestion 
institucional en la alcaldfa menor. Para su deterrninacion se definieron los siguientes 
componentes: 

2.3.1 Camas hospital arias, XR,: Definido como el numero de camas normalizado del conjunto de 
hospitales e instituciones de salud en el area que cubre la alcaldfa menor. 

2.3.2 Recurso humano en salud , Xr 2 : Definido como el numero normalizado de medicos y 
enfermeras que trabajan en las instituciones de salud local izadas en la jurisdiccion de la 
alcaldfa menor. 

2.3.3 Espacio publico , Xr 3 : Definido como el area normalizada de espacio disponible para el 
albergue o alojamiento temporal y la atencion masiva de emergencias en el area que cubre la 
alcaldfa menor. 

2.3.4 Personal de rescale , Xr 4 : Definido como el numero normalizado de socorristas de Cruz Roja, 

Defensa Civil y Bomberos disponible en la zona de influencia de la alcaldfa menor. 

2.3.5 Nivel de desarrollo , Xr 5 : Definido como la calificacion del nivel de calidad de vida. 
organizacion y planeamiento urbano en el area de la alcaldfa menor, valorado por la 
Secretana de Planeacion de la ciudad. 
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2.3.6 Operatividaci en emergencies , Xr 6 : Definida como la calificacion de la preparacion y la 
capacidad de respuesta institucional del comite operativo de emergencias de la alcaldfa 
menor, valorado la Oficina de Prevencion y Atencion de Emergencias de la ciudad. 

La figura 4 presenta los valores de los indicadores de exposicion, fragilidad social y falta de 
resiliencia que componen el descriptor de vulnerabilidad del contexto. La figura 5 ilustra los valores 
agregados del mismo descriptor. 




Figura 4 Descriptor de vulnerabilidad del contexto desagregado en indicadores. 
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Figura 5 Descriptor de vulnerabilidad del contexto. 






Las figura 6 presenta los valores de los descriptores de amenaza y vulnerabilidad del contexto, la 
figura 7 ilustra los valores de los indices de rieSgo sismico ffsico y del contexto y la figura 8 
presenta los valores del fndice de riesgo sismico total ordenados de menor a mayor. Este tipo de 
graficas permiten categorizar y priorizar las alcaldfas que presentan los mayores valores 
comparativos. Por otra parte, el conjunto de graficas de la figura 9 presenta la desagregacion de los 
indicadores del riesgo sismico ffsico, de descriptor de amenaza sfsmica y de los indicadores de 
exposicion, fragilidad social y falta de resiliencia de una de las alcaldfas menores. Este tipo de 
graficos, entre otros, permite identificar que indicadores tienen mayor incidencia comparativa y por 
lo tanto a cuales de ellos se deben dirigir las medidas de mitigacion y prevencion. 
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Figura 6 Descriptores de amenaza y vulnerabilidad del contexto. 




Figura 7 Indices de riesgo sfsmico ffsico y riesgo sfsmico del contexto. 




Figura 8 Categorizacion de las alcaldi'as menores por el l'ndice de riesgo sfsmico total. 
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Escenario de perdidas por terremoto. 





Amenaza sfsmica 



Exposicion 




Fragilidad social 



Resiliencia 



Figura 9 Desagregacion de indicadores para una alcaldia menor. 
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3 CONCLUSION 

Se ha desarrollado un modelo para la estimacion de un indice de riesgo sismico relative de las 
localidades que conforman un centra urbano metropolitano. Esta metodologfa, aplicada a Santa Fe 
de Bogota Colombia, ha permitido categorizar la amenaza, la vulnerabilidad y el riesgo sismico de 
las diferentes alcaldias menores de la ciudad desde una perspectiva holfstica. En este metodo se 
puede actualizar facilmente el valor de las variables, lo que favorece la realizacion de analisis de 
sensibilidad y calibracion. De la misma manera, puede llevarse a cabo el seguimiento del eScenario 
de riesgo y de la efectividad y eftciencia de las medidas de prevencion y mitigacion. Una vez 
graftcados los resultados, para cada una de las localidades, es sencillo identificar los aspectos mas 
relevantes del riesgo sismico relativo, sin la necesidad de realizar mayores esfuerzos de analisis e 
inteipretacion de resultados. La principal ventaja de la tecnica aqui planteada es la posibilidad de 
“devolverse” mediante la desagregacion de los indices en descriptores y estos a su vez en 
indicadores. e identificar porque razon una localidad presenta un mayor indice de riesgo. Esta virtud 
del metodo permite verificar los resultados y priorizar las acciones prevencion y planificacion que 
se.deben implementar para intervencion y modificacion de las condiciones que mas influyen en el 
riesgo sismico de la ciudad. 
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APENDICE A. 

RIESGO SISMICO FISICO DE SANTA FE DE BOGOTA 

Este apendice presenta un resumen de la base conceptual, los analisis, los resultados y las 
aplicaciones de los estudios de microzonificacion sfsmica y de la estimacion de los escenarios de 
perdidas en caso de sismos probables, establecidos por los estudios de amenaza sfsmica y el 
comportamiento de los suelos en el area donde se encuentra localizada la ciudad. 

MICROZONIFICACION SISMICA 

La ciudad de Bogota, Distrito Capital de Colombia, esta localizada en las inmediaciones de un 
ambiente sismotectonico que demuestra, desde el punto de vista geologico, actividad sfsmica 
reciente, aunque desde el punto de vista historico no se tenga la percepcion de dicha actividad 
debido a la baja frecuencia relativa de eventos importantes. Esta aparente contradiccion debe verse 
con extrema cautela dado que no es correcto suponer que como en los ultimos arios no ha ocurrido 
un sismo intenso asf han de seguir las cosas. El sismo del 19 de enero de 1995, de magnitud 6.5 
grados Richter que ocurrio a 120 km de Bogota, tuvo una aceleracion de solo el 0.03g (despues 
de amplilicado por el suelo), es decir 7 veces menos que el valor maximo probable esperado 
(0.20% g) para el cual se deben disenar los edificios en la ciudad; no obstante este evento causo 5 
vfctimas, 28 heridos, dafios en 40 edificios, crisis en los sistemas de comunicaciones, 30 
accidentes de transito y dos incendios en la capital. 

No obstante que existfa un estudio general de amenaza sfsmica de Colombia realizado en 1983. 
el cual fue recientemente actualizado (Asociacion Colombiana de Ingenierfa Sfsmica, 1996), para 
este estudio se llevo a cabo un nuevo analisis detallado de la amenaza sfsmica regional y local de 
Bogota. El Ingeominas y la Universidad de los Andes realizaron un trabajo conjunto de sismologfa 
historica, tectonica regional y neotectonica con el fin de precisar las trazas y la actividad de las 
fallas geologicas mas importantes del area en donde se encuentra la ciudad. Igualmente, se llevo a 
cabo una revision y ajuste del catalogo sismico de eventos, el cual ha sido alimentado de manera 
notable en los ultimos anos con sismos registrados tanto en acelerometros digitales como en las 
estaciones sismologicas de la nueva red sfsmologica nacional que opera via satelite desde 1991 . 

Del estudio de amenaza sfsmica se pudo obtener que Bogota se encuentra amenazada 
principalmente por el conjunto de fallas que conforman la falla Frontal de la Cordillera Oriental, las 
cuales pueden llegar a tener eventos del orden de 7 a 7.5 grados en la escala de magnitudes de 
Richter a una distancia del orden de 60 km al oriente de la ciudad. Sin embargo, la existencia de 
otras fallas activas alrededor y el potencial de ocurrencia de sismos lejanos que pueden amplificarse 
notablemente debido a las propiedades del suelo de la ciudad, ratifica que la amenaza sfsmica de 
Bogota no depende de una sola fuente sismogenica sino de una amplia provincia de accidentes 
geotectonicos. 

Teniendo en cuenta la informacion existente y complementandola con nuevos estudios 
geologicos, hidrogeologicos, geotecnicos y geoffsicos (gravimetria, microtremores, down hole, 
refraction, reflexion), se propuso un mapa de microzonificacion sismica, tanto de la zona plana de 
origen lacustre como de los cerros orientales de la ciudad, de acuerdo con los modelos de respuesta 
sfsmica unidimensional y bidimensional evaluados y calibrados con el programa SHAKE 91 y 
mediante modelos de elementos finitos utilizando el programa generador de mallas ANSYS y 
QUAD4M (Idriss et al, 1992). Esta modelacion se obtuvo utilizando registros de aceleracion reales 
y artificiales en roca y en el suelo blando de la ciudad. Se realizaron cerca de 117 mediciones de 
microtremores, lo cual permitio obtener un mapa de perfodos de vibration predominante que se 
ilustra en la figura 2 (Yamfn, Ojeda, 1995) y 38 perforaciones profundas entre los 20 y los 200 m 
que permitieron caracterizar, mediante ensayos dinamicos de laboratorio, los tipos de suelos de 
Bogota. Se encontro que estos suelos, que tienen un comportamiento elastico sorprendente, 
amplifican entre 3 a 10 veces, o mas, la aceleracion maxima registrada en la roca. Igualmente, se 
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encontro un efecto de amplification por topografi'a notable que puede sigmficar la ocurrencia de 
deslizamientos en algunos sitios en los cerros orientales de la ciudad. 

La figura A. 1 corresponde a la microzonificacion sismica, en la que se diferencian cinco zonas 
para las cuales se determinaron los espectros de disefio con los cuales se deben analizar y disefiar las 
edificaciones en cada zona, figura A.2. Utilizando esta metodologia se obtuvieron .gualmente los 
valores de aceleracion espectral para otros sismos de menor intensidad, con el hn de estimar otros 
escenarios de perdidas potentiates en la ciudad en caso de sismos moderados. La metodologia 
desarrollada para la determination de los danos en edificaciones y lrneas Vitales parte del 
conocimiento de la respuesta espectral para cada escenario si'smico hipotetico, en cada una de las 
zonas identificadas en la microzonificacion sismica de la ciudad. 



VULNERABILIDAD FISICA URBANA 

El analisis del potential de danos se realizo llevando a cabo una identilicacion del tipo de 
edificaciones que se han construido en la ciudad desde su fundacion. Para el efecto fue necesario 
hacer un a revision del crecimiento historico de Bogota y de la manera como las tecnologi'as de la 
construction se fueron modificando con el pasar del tiempo. Por otra parte y teniendo en cuenta que 
las tipologi'as de construction tambien cambian de acuerdo con la capacidad economica, fue 
necesario hacer un analisis de los estratos socioeconomicos y de la informacion que en relation con 
las zonas de uso, limitation de altura de las edificaciones y zonas de tratamiento tiene la Secretaria 
de Planeacion de la ciudad. Finalmente, con base en informacion sistematizada del Catastro Distntal 
y una amplia verification en el campo realizada por cuadrillas de evaluadores se logro desarrollar 
una base de information que permitiera conocer que tipo de tipologi'as de edificaciones son las mas 
comunes en la ciudad y la manera como estan distribuidas en cada una de las manzanas, barrios y 
sectores de la ciudad. 

Teniendo en cuenta que las propiedades dinamicas y las caracteristicas de resistencia, rigidez y 
disipacion de energi'a inelastica de las edificaciones ante los terremotos vari'an notablemente 
dependiendo del tipo de material del sistema estructural, de la altura, de la tecnica de construccion, 
de la edad, y de otros factores, se consideraron las tipologi'as indicadas en la Tabla 1 . 



Tabla A.l Tipos de edificaciones. 



TIPOS DF. F.STRUCTURAS 

1 . Mamposteria, 1-2 pisos 

2. Mamposteria >= 3 pisos 

3. Porticos de CR hasta 5 pisos 

4. P6rticos de CR 5- 10 pisos 

5. Porticos de CR mas de 10 

6. Industrial liviana 

7. Industrial pesada 

8. Informal liviana 

9. Informal pesada 

10. Otras 



USO 

1 . Residencial 

2. Comercial 

3. Education 

4. Salud 

5. Institucional 

6. Industrial 

7. Otro 



CUBIERTA 

1. Liviana 

2. Teja de barro 

3. Placa 



ANQ DE CONSTRUCCION 

1 . Post codigo 

2. Moderna 

3. Premodcrna 

ENTREPISOS 

1 . Concreto 

2. Madera 

3. Otro 

MANTENIMIENTO 

1. Bueno 

2. Regular 

3. Malo 



Es importante notar que este tipo de clasificacion se definio con base en las caracteristicas 
arquitectonicas y constructivas observadas del conjunto general de edificaciones de la ciudad. No 
menos del 85% de la poblacion se encuentra alojada en edificaciones de uno y dos pisos tfpicamente 
de mamposteria. Un amplio numero de los conjuntos multifamiliares tienen como li'mite superior los 
cinco pisos debido a la reglamentacion urbana y a los requerimientos de ascensores en los edificios. 
os edificios con estructuras metalicas no se consideran como una tipologi'a adicional, ya 
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Figura A1 Mapa de Microzonificacion Sismica de Santa Fe de Bogota 
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ESPECTRO DE DISEftO Y RESPUESTA. 
ZONA1 - CERROS 
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Figura A2 Espectros de aceleracion para las diferentes zonas 
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que al momento de llevar a cabo este trabajo solo se habfa construido en la ciudad cerca de 30 
edificios de porticos en acero. 

El levantamiento de la informacion consistio en identificar en cada una de las manzanas de la 
ciudad las tipologfas de las edificaciones, determinar la cantidad de area construida aproximada de 
las mismas y sus principales caracterfsticas. Utilizando los pianos mas recientes de la ciudad, 
elaborados por el Instituto Geografico Agustfn Codazzi a escala 1:2000, y con la ayuda de una 
planilla debidamente disenada para el efecto se identificaron, la manzana, el numero de la 
edificacion y sus caracterfsticas. Para efectos de determinar la altura de las edilicaciones se 
utilizaron fotograffas aereas laterales similares a las publicadas en el libro "Bogota desde el Aire", 
(Horney J. et al , 1994), las cuales permitfan conocer la altura de las edificaciones por el angulo en 
que fueron tomadas, a diferencia de las fotograffas aereas tfpicas. Igualmente se utilizo informacion 
de la Secretarfa de Planeacion y el Catastro Distrital, y finalmente de las visitas de los grupos 
tecnicos de verificacion que recorrieron la ciudad, realizando inspecciones manzana por manzana en 
las zonas donde la heterogeneidad de las tipologfas asf lo ameritaban. 

En total se cubrio una area de 405 km 2 , correspondiente a 46.000 manzanas sobre un perfmetro 
de la ciudad del orden de 92 km. El area construida total sumo 317 millones de metros cuadrados, 
donde las estructuras informales y de mamposterfa son el 66%, los porticos de concreto retorzado 
hasta de 5 pisos son el 20% y los porticos mayores a 5 pisos son el 9%. Una vez obtenida la 
informacion por manzanas se llevo a cabo una sectorizacion de la ciudad por "celdas", las cuales 
pasaron a ser las unidades de analisis del estudio. Esto con el fin de presentar la informacion 
adecuadamente mediante un S1G. En cada celda la aceleracion del suelo se considero constante y se 
trato de mantener cierta homogeneidad de acuerdo con la distribucion de las tipologfas estructurales; 
ademas, sus perfmetros se hicieron coincidir en su mayor parte con lfmites ffsicos existentes tales 
como vfas principales, rios y parques, como tambien con los lfmites de la sectorizacion de las 
alcaldfas menores de la ciudad. 

Con el fin de determinar los diferentes niveles de dano que puede presentar una edificacion, se 
adopto la metodologfa ya ampliamente generalizada propuesta en el ATC-13 (Applied Technology 
Council, 1985) basada en Estados de Dano.Estos valores de probabilidad de dano para cada valor de 
intensidad sfsmica propuestos por el metodo ATC-13, son basados en las matrices de probabilidad 
de danos desarrolladas por Whitman et al (1973), aplicando la tecnica estadfstica del Metodo 
Delphi, en el cual se utilizo el criterio y las apreciaciones de un amplio grupo de expertos. Dichas 
matrices de probabilidad permiten deducir unas funciones de vulnerabilidad para ciertos valores de 
intensidad, en las cuales cada estado de dano puede expresarse en terminos ffsicos. 

De acuerdo con el criterio y experiencia obtenida por expertos locales de sismos colombianos y 
utilizando la ecuacion la Trifunac y Brady (1975) que relaciona la Intensidad MM y la aceleracion 
pico del suelo y llevando a cabo una revision de las diferentes funciones de vulnerabilidad 
desarrolladas por varios autores como Algermissen (1984) y Steinbrugge (1982), Sauter y Shah 
(1978), Petrovsky y Milutinovic (1985), Ordaz et al.(1994) se ajustaron las matrices de danos 
propuestas por el ATC-13 para el caso de las tipologfas de edificacion definidas para Bogota. 

El valor de la amenaza sfsmica en cada zona de la ciudad esta representado por un espectro de 
aceleracion obtenido para cada sismo postulado. En consecuencia, para conocer el valor de 
aceleracion espectral al cual esta sometida cada tipologfa de edificacion fue necesario estimar el 
periodo fundamental de vibracion de cada una. Ahora bien, utilizando el metodo desarrollado por 
Scholl et al (1982) para relacionar valores de aceleracion espectral e intensidad teniendo en cuenta 
tres tipos de suelo similares a los que se encuentran en Bogota y 546 acelerogramas de terremotos 
registrados en todo el mundo desde 1933 a 1979, se realizo un procedimiento de correlacion que 
permitiera mediante una funcion de transferencia, conocidos los valores de aceleracion espectral, 
utilizar las funciones de vulnerabilidad modificadas, que fueron desarrolladas para valores de 
intensidad en la escala de Mercalli Modificada. El costo de reposition se obtuvo para cada tipologfa 
estructural con base en el costo por metro cuadrado de construction dependiendo del tipo de estrato 
economico. Dicho valor vario desde US$ 240 hasta US$ 864 dolares promedio por metro cuadrado, 
para edificaciones en mamposterfa en estrato bajo y porticos de concreto reforzado de mas de diez 
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pisos en estrato alto respectivamente. De esta informacion se obtienen las perdidas economicas 
directas como un porcentaje del costo de reposition para cada tipologi'a de edification. Debido a la 
falta de informacion local, para efectos de estimar la cantidad de personas que pueden presentar 
heridas e incluso la muerte, se utilizaron estimativos similares a los propuestos por Whitman et al. 
(1973) y ajustados en el ATC-13, que son en general aceptados para los pai'ses oecidentales. Las 
figuras A. 3 y A.4 ilustran el area de construction danada considerando los dos escenarios mas 
graves de amenaza sismica para la ciudad. 

Teniendo en cuenta la informacion disponible aportada por las diferentes empresas de servicios 
publicos involucradas, se desarrollo un modelo simplificado que permitiera relacionar de mantra 
global, la probabilidad de fallo o dafio de los componentes principales de cada sistema con la 
intensidad de los terremotos probables en la ciudad. Para el estudio de las principales redes de 
distribution se utilizo la metodolpgi'a de estimation de dafio propuesta por el ATC-13. El dafio se 
evaluo como un porcentaje de la longitud cafda de cable para el caso de lfneas de alia tension y poi 
numero de roturas por km para el caso de las tuberias de acueducto y gas. Para el calculo de la 
vulnerabilidad de las edificaciones principales de los sistemas de energi'a electrica (subestaciones 
electricas) y telefonos (centrales telefonicas) se utilizo la metodologi'a aproximada del ATC-21 . 

ESCENARIOS DE PERDIDAS 

Se consideraron tres escenarios de perdidas: un sismo lejano con una aceleracion en el 
basamento rocoso de 0.038g, un sismo cercano moderado de aceleracion de 0.1 2g y un sismo 
cercano fuerte de aceleracion de 0.20g en el basamento rocoso. Del trabajo realizado se puede 
concluir que aun cuando la amenaza sismica para de Bogota no corresponde a valores extremos en 
el pals, el riesgo si lo es debido al alto grado de vulnerabilidad de sus edificaciones, las cuales hasta 
hace muy pocos aiios fueron construidas sin tener en cuenta criterios sismorresistentes. Solo a partir 
de 1984, afio en el cual se aplico el primer estudio de amenaza sismica en el pais para efectos de 
expedir por primera vez una normativa sismica obligatoria, se podrfa decir que las edilicaciones 
cuentan con algiin grado de protection en caso de terremoto; obviamente siempre que hayan sido 
construidas en forma correcta. Infortunadamente, con anterioridad a 1984 no se tuvieron 
requerimientos para soportar solicitaciones dinamicas y cargas laterales, razon por la cual las 
edificaciones disenadas y construidas antes de ese afio en su mayoria son altamente vulnerables, 
como lo han demostrado sismos incluso mucho menores a los de disefio en otras zonas del pais. De 
acuerdo con la informacion obtenida de las evaluaciones presentadas anteriormente, se puede 
resumir en algunas cifras el alto riesgo sismico de Bogota, entendido como el potencial de 
consecuencias sociales, ambientales y economicas en caso de terremoto. La Tabla A..2 presenta las 
principales cifras obtenidas del estudio para el inventario de las edificaciones y para cada uno de los 
escenarios sismicos hipoteticos considerados. 

Tabla A.2 Resultados de Perdidas Totales - Escenarios Sismicos 



ESCENARIO 

SISMICO 

HIPOTET1CO 


EN EL DIA 


EN LA NOCHE 


SIN VIVIENDA 
(10% de Afectados 
Habitantes) 


AREA 

DESTRUIDA 
(Millones m 2 ) 


COSTO 

(US$ 

Millones) 


Muertos 


Heridos 


Muertos 


Heridos 


SISMO CERCANO 
FUERTE (0.20&) 


3500 


20000 


4500 


26000 


74000 


33.8 


•14000 


SISMO CERCANO 

MODERADO 

(0.12s) 


1600 


9000 


1400 


7700 


44000 


20.6 


8800 


SISMO LEJANO 
(0.03s) 


300 


1600 


350 


1900 


27000 


12.3 


5100 
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Figura A.3 Porcentaje de Area destruida para el sismo Cercano Moderado 




Figura A.4 Porcentaje de Area destruida para el sismo Cercano Fuerte 
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Por dificultades de acceso e interrupcion del transito, lo cual impedirfa el rescate y la atencion 
medica inmediata, se sumaran los heridos muy graves a los valores estimados de muertos obtenidos 
de la Tabla A.2. Igualmente, teniendo en cuenta que un estudio de esta naturaleza no puede ser 
exacto, se aproximaran las cifras obtenidas del calculo a numeros globales. que sirven como valores 

indicatives de la dimension de la catastrofe. 

En la eventual idad de la ocurrencia de un sismo en la Falla Frontal de la Cordillera Oriental, que 
cause una aceleracion del orden de 0.1 2g a nivel del basamento rocoso, el cual se considera podna 
ser un sismo con un perfodo de retomo de 100 anos, similar a los sismos que afectaron a la c.udad 
en 1785, 1827 y 1917, se esperaria una fuerte respuesta en la zona oriental de la ciudad cercana a los 
cerros, lo que afectarfa notablemente edificaciones localizadas en las zonas de suelos rocosos y de 
piedemonte. Este sismo, a diferencia del anterior, causaria graves danos estructurales y no 
estructurales en edificaciones de diferentes alturas. Se estima que si el sismo ocurre durante el dm 
habria del orden de 1.600 muertos y 9.000 heridos. Si ocurre en la noche habria cerca de 1.400 
muertos y 7 700 heridos. Habria danos equivalentes a la destruccion del 6.8% del area construida de 
la ciudad, lo que corresponde a 20.6 millones de metros cuadrados de edificaciones. cuyos costos 
podrian alcanzar los 8.800 millones de dolares. Se estima que en el area afectada habria 435.000 
habitantes, de los cuales se estima que no menos del 10% podrian tener problemas de alojamiento, 
es decir, del orden de 44.000 personas sin vivienda. Una fuerte respuesta en la zona oriental de la 
ciudad cercana a los cerros, afectarfa notablemente segmentos de la red del acueducto localizados en 
las zonas de suelos rocosos y de piedemonte. No obstante, en la zona plana habria danos 
importantes en todos los sistemas de redes. La red del acueducto se estima podrfa tener en total, en 
los diferentes tipos de tuberfa, entre 300 y 350 rupturas, la red de gas natural entre 25 y 30 danos 
importantes y la red de alta tension del orden de 4 kilometres y medio de lfneas de cables caidos. 
Varias subestaciones de energfa presentarfan danos que contribuin'an a interrumpir los servicios. 
Esto agravado por el amplio numero de vfetimas, los posibles incendios por los escapes adicionales 
de gas propano, significant una grave crisis para la ciudad y la nacion, pues el tiempo de 
recuperacion de los servicios y la atencion a la poblacion serfa notablemente deficiente, razon por 
la cual las medidas de reduction de vulnerabilidad y riesgo de los sistemas como la elaboracion de 
planes de contingencia acorde con este escenario son acciones de especial tmportancia, mas si se 
tiene en cuenta la alta probabilidad de ocurrencia de un evento como este, que ya ha ocurndo en el 

pasado en la ciudad. , , _ OA , 

Ante un sismo que se presente en la misma falla que genere aceleraciones del orden de 0.20gjo 
cual equivale a un sismo de 475 anos de periodo de retomo, considerado como el sismo de diseno 
para el cual se deben diseiiar las estructuras de las edificaciones que se construyen en la ciudad, se 
tendrfa el escenario de perdidas mas notable de los considerados en este estudio. Las perdidas de 
nuevo se presentarfan de manera intensa en las cercanfas de los cerros orientales, pero habria una 
participation elevada de otros sectores al norte y al sur. Este sismo causaria graves danos 
estructurales y no estructurales en todos los tipos de edificaciones. Habria colapsos totales de 
estructuras en diferentes sitios de la ciudad, con menor incidencia en la zona occidental donde los 
danos en general serian menores. Se estima que si el sismo ocurre durante el dfa habria del orden de 
3.500 muertos y 20.000 heridos. Si ocurre en la noche habria 4.500 muertos y 26.000 heridos. 
Habria danos equivalentes a la destruccion del 10.9% del area construida de la ciudad. lo que 
corresponde a 33.8 millones de metros cuadrados de edificaciones. cuyos costos serian de 14.000 
millones de dolares. Se estima que en el area afectada habria 738.000 habitantes, de los cuales se 
estima que no menos del 10% podrian tener problemas de alojamiento, es decir. del orden de 74.000 
personas sin vivienda. La red del acueducto se estima podrfa tener entre 450 y 500 rompimienios, la 
red de gas natural del orden de 60 danos importantes y la red de alta tension cerca de 6 kilometres v 
medio de lfneas de cables colapsadas. Un amplio numero de subestaciones de energfa del norte, el 
sur y el borde oriental, y la mayorfa de las centrales telefonicas del norte el centro-oriente y sur de la 
ciudad presentarfan danos importantes para el funcionamiento de los servicios. Serfa la mayoi crisis 
factible para la nacion debido a la concentration de poblacion, bienes y servicios en la capital y por 
lo tanto no solo se debe hacer un esfuerzo notable de las instituciones para identificar recursos 
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financiers y presupuestales cada ano para realizar medidas efectivas de reduccion de 
vulnerabilidad y riesgo de los sistemas, sino la elaboracion de un plan especffico de contingencia 
para responder de la mejor manera posible ante este escenario extremo, no solo considerando la 
atencion a la poblacion inmediatamente y semanas despues del evento sino para efectos de 
rehabilitacion y recuperacion de la infraestructura afectada. 

Estas cifras no son exactas y solamente permiten dimensionar en ordenes de magnitud la 
problematica que para Bogota y Colombia habrfa en la eventualidad de la ocurrencia de cualquiera 
de los escenarios hipoteticos propuestos, que obviamente tienen una probabilidad de ocurrencia y 
pueden Uegar a ocurrir segun la informacion disponible y las evidencias existentes en la actualidad. 
Por esta razon, es fundamental incorporar este tipo de evaluaciones a los programas de 
ordenamiento urbano y definir procedi mientos y protocolos de respuesta institucional. areas de 
alojamiento temporal y ejercicios de simulacion de eventos hipoteticos. 

Este apendice es un resumen del artfculo realizado por Cardona y Yamfn (1997), el cual describe 
los estudios realizados por la Universidad de los Andes y el lngeominas, entidades que hicieron un 
especial esfuerzo tecnico-cientffico y economico para la estimacion del riesgo sfsmico de Santa Fe 
de Bogota, el cual conto con el apoyo financier e institucional de la Direccion Nacional para la 
Prevencion y Atencion de Desastres y a la Unidad Distrital de Prevencion y Atencion de 
Emergencias de Bogota. 
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RESUMEN 

A partir de los trabajos de Hanks y McGuire (1981), Joyner y Boore (1984), Luco (1985) y 
Ordaz (1992) se ha desarrollado procedimientos para estimar las aceleraciones maximas del suelo 
con base en modelos sismologicos del espectro radiado. En el presente trabajo, esto se repite para 
diversas distancias y magnitudes, y las aceleraciones calculadas se comparan con las registradas 
durante sismos colombianos, con el fin de determinar los parametros relevantes del modelo 
sismologico del espectro y asi poder llegar leyes de atenuacion de aceleracion maxima del suelo que 
cumplen con los requisitos de desviacion y error medio mfnimos. Se hace ver que el inetodo 
presentado es poderoso, ya que con el pueden tambien derivarse leyes de atenuacion para otras 
cantidades de interes para la ingenieria sfsmica, como velocidad, desplazamiento u ordenadas 
espectrales, mediante simples relaciones teoricas. Asimismo, mediante el uso de funciones de 
transferencia, empfricas o teoricas, se pueden generar leyes de atenuacion para sitios especificos, 
incluyendo los efectos locales. Se hace una comparacion con el reciente sismo del Quindio, 
Colombia, para estaciones situadas sobre terreno firme y se comprueba el buen comportamiento de 
la ley propuesta. Ademas, se encuentran buenas aproximaciones para ordenadas espectrales. Se 
muestra que este enfoque es especialmente util para regiones del mundo en que no se tienen bases 
de datos amplias de registros de movimiento fuerte. 



ABSTRACT 

Starting with the works of Hanks and McGuire (1981), Joyner and Boore (1984), Luco (1985) 
and Ordaz (1992), procedures have been developed to estimate peak ground accelerations from 
seismological models of the radiated spectra. In this work, this procedure is applied for several 
distances and magnitudes. Computed accelerations are compared with recorded values in Colombia, 



©1999 ESPE, Quito, Ecuador 



ISSN: 1390-0315 



Recibido: Marzo de 1999 



